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CN ■ Abstract 

\ In this paper, a class of relay networks is considered. We assume that, at a node, outgoing channels to its 

■ neighbors are orthogonal, while incoming signals from neighbors can interfere with each other. We are interested in 

the multicast capacity of these networks. As a subclass, we first focus on Gaussian relay networks with interference 
and find an achievable rate using a lattice coding scheme. It is shown that there is a constant gap between our 
achievable rate and the information theoretic cut-set bound. This is similar to the recent result by Avestimehr, 
Diggavi, and Tse, who showed such an approximate characterization of the capacity of general Gaussian relay 
networks. However, our achievability uses a structured code instead of a random one. Using the same idea used in 
the Gaussian case, we also consider linear finite-field symmetric networks with interference and characterize the 
Q I capacity using a linear coding scheme. 
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■ I. Introduction 

O \ Characterizing the capacity of general relay networks has been of great interest for many years. In this 
^ paper, we confine our interest to the capacity of single source multicast relay networks, which is still 
^ an open problem. For instance, the capacity of single relay channels is still unknown except for some 
• ^ special cases [IJ. However, if we confine the class of networks further, there are several cases in which 
^ the capacity is characterized. 

c3 ■ Recently, in [[2l], the multicast capacity of wireline networks was characterized. The capacity is given by 
the max-flow min-cut bound, and the key ingredient to achieve the bound is a new coding technique called 
network coding. Starting from this seminal work, many efforts have been made to incorporate wireless 
effects in the network model, such as broadcast, interference, and noise. In [3J, the broadcast nature was 
incorporated into the network model by requiring each relay node to send the same signal on all outgoing 
channels, and the unicast capacity was determined. However, the model assumed that the network is 
deterministic (noiseless) and has no interference in reception at each node. In [4J, the work was extended 
to multicast capacity. In [5J, the interference nature was also incorporated, and an achievable multicast rate 
was computed. This achievable rate has a cut-set-like representation and meets the information theoretic 
cut-set bound Il27ll in some special cases. To incorporate the noise, erasure networks with broadcast or 
interference only were considered in Q, [[8]]. However, the network models in [7J, [8 J assumed that the 
side information on the location of all erasures in the network is provided to destination nodes. Noisy 
networks without side information at destination nodes were considered in LI 2.1 and LI 3,1 for finite-field 
additive noise and erasure cases, respectively. 
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Along the same lines of the previous work on wireless networks mentioned above, we consider the 
multicast problem in a special class of networks called relay networks with interference. More specifically, 
we assume that all outgoing channels at each node are orthogonal, e.g., using frequency or time division 
multiplexing, but signals incoming from multiple neighbor nodes to a node can interfere with each other. 
Since wireless networks are often interference limited, our setup focuses on the more important aspect 
of them. This model covers those networks considered in 0], lITOll , 111211 . Our interest in the relay 
networks with interference was inspired by [14J, in which the capacity of single relay channels with 
interference was established. In this paper, we focus on two special subclasses of general networks with 
interference; Gaussian relay networks with interference and linear finite-field symmetric networks with 
interference. 

For the Gaussian relay networks with interference, we propose a scheme based on nested lattice codes 
lfT9ll which are formed from a lattice chain and compute an achievable multicast rate. The basic idea of 
using lattice codes is to exploit the structural gain of computation coding [11 J, which corresponds to a 
kind of combined channel and network coding. Previously, lattices were used in Gaussian networks in 
ifTOll . and an achievability was shown. However, our network model differs from the one in HI Oil in that 
we assume general unequal power constraints for all incoming signals at each node, while an equal power 
constraint was mainly considered in [lOJ. In addition, our lattice scheme is different from that in [lOJ in 
that we use lattices to produce nested lattice codes, while lattices were used as a source code in [lOJ. 

We also show that our achievable rate is within a constant number of bits from the information theoretic 
cut-set bound of the network. This constant depends only on the network topology and not on other 
parameters, e.g., transmit powers and noise variances. This is similar to the recent result in [6J, which 
showed an approximate capacity characterization for general Gaussian relay networks using a random 
coding scheme. However, our achievability uses a structured code instead of a random one. Thus, our 
scheme has a practical interest because structured codes may reduce the complexity of encoding and 
decoding. 

Finally, we introduce a model of linear finite-field symmetric networks with interference, which gen- 
eralizes those in lfT2ll . ifTSll . In the finite-field case, we use a linear coding scheme, which corresponds to 
the finite-field counterpart of the lattice coding scheme. The techniques for deriving an achievable rate 
for the finite-field network are basically the same as those for the Gaussian case. However, in this case, 
the achievable rate always meets the information theoretic cut-set bound, and, thus, the capacity is fully 
established. 

This paper is organized as follows. Section [III defines notations and parameters used in this paper and 
introduces the network model and the problem of interest. In Section [IIIl we analyze Gaussian relay 
networks with interference and give the upper and lower bounds for the multicast capacity. In Section [IVl 
we define a model of linear finite-field symmetric networks with interference and present the multicast 
capacity. Section [V] concludes the paper. 

II. Relay networks with interference 
A. System model and notations 

We begin with a description of the class of networks that will be dealt with in this paper. The memoryless 
relay networks with interference are characterized such that all outgoing channels from a node to its 
neighbors are orthogonal to each other. We still assume that incoming signals at a node can interfere with 
each other through a memoryless multiple-access channel (MAC). An example of this class of networks is 
shown in Fig. [H Some special cases and subclasses of these networks have been studied in many previous 
works W\. f9l, fT05, fT3l, [fTl. 

We will begin by giving a detailed description of the network and some definitions of the parameters. 
The network is represented by a directed graph Q = {V,E), where l/ = {l,...,|y|}isa vertex set and 

C y X y is an edge set. Each vertex and edge correspond to a communication node and a channel in 
the network, respectively. In this paper, we focus on a multicast network: vertex 1 represents the source 
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node and is denoted by s, and the set of destination nodes is denoted by D, where s ^ D. It will be 
assumed that the source node has no incoming edge, and the destination nodes have no outgoing edge. 
All the other nodes, which are neither the source nor the destination, are called the relay nodes. Since all 
broadcast channels in the network are orthogonal, we associate a discrete or continuous random variable 
Xul at time t with edge {u, ^) G as a channel input (output of a node). As a channel output (input 
of a node), we associate a discrete or continuous random variable yj^^ at time t with node v eV \ 
From now on, we sometimes drop the superscript when doing so causes no confusion. 
At node v eV, tht set of incoming and outgoing nodes are denoted by 

A(^) = {u : {u,v) G E}, 
Q{v) {w : {v,w) G E} . 

Set S C V is called a cut if it contains node s and its complement contains at least one destination 
node d e D, i.e., S^HD ^ Let F denote the set of all cuts. The boundaries of S and are defined as 

S = {u :3v s.t. {u, v) e E,u e S,v e S""} , 
S"" = {v : 3u s.t. {u, v) e E,u e S,v e S""} . 

For node v G S'^, the set of incoming nodes across S is defined as 

As{v) = A{v) n S = A{v) n S, 
For any sets Si CV and S2 ^V, wt define 

-^^1,^2 = {^u,v ' {u, v) e E,u e Si,v e S2} , 

and 

Xa{v) = {Xu,v ' u G A(^)}. 

Using the aforementioned notations, we can formally define the class of networks of interest. The 
memoryless relay network with interference is characterized by the channel distribution function 

p {yv\xvy) = p (2/2|^A(2)) P {ys\xA(3)) '"P (2/m|^a(m)) 
over all input and output alphabets. 
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B, Coding for the relay network with interference 

The multicast over the relay network consists of encoding functions ful{-), {u,v) e E, t = 1, . . . , N, 
and decoding functions gd{'), d G D. The source node s has a random message W that is uniform over 
{1, . . . , M} and transmits 

X(t) = fit) (w) 

at time t on the outgoing channels {s,w), w G Q{s). The relay node v transmits 

X(t) = fit) (Y^-^) 

at time t on the outgoing channels {v,w), w G B(^), where Y^~^ = (Yv^\ . . . ,1"^^"^^^ At destination 
node d e D, after time A^, an estimate of the source message is computed as 

Then, the probability of error is 

Pe = Pr|^U^{^.(r/)^l^}}. (1) 

We say that the multicast rate R is achievable if, for any e > and for all sufficiently large A^, encoders 
and decoders with M > 2^^ exist such that Pg < e. The multicast capacity is the supremum of the 
achievable multicast rates. 

As stated in Section [D we are interested in characterizing the multicast capacity of the memoryless 
relay networks with interference. However, as shown in [13J, even for a relatively simple parallel relay 
channel, finding the capacity is not easy. Thus, we further restrict our interest to the Gaussian networks 
in Section [Illl and the linear finite-field symmetric networks in Section [IVl 

III. Gaussian relay networks with interference 

In this section, we consider Gaussian relay networks with interference. At node v at time t, the received 
signal is given by 

y W = \^ I zi^) 

ueA{v) 

where Zy'' is an independent identically distributed (i.i.d.) Gaussian random variable with zero mean and 
unit variance. For each block of channel input (^X^l, . . . , X^^i^, we have the average power constraint 
given by 

t=l 

In fTO'], Nazer et al. studied the achievable rate of the Gaussian relay networks with interference for the 
equal power constraint case, where Pu^y = Py for all u e A{v). In our work, we generalize it such that 
Pu^y's can be different. The main result of this section is as follows. 

Theorem 1: For a Gaussian relay network with interference, an upper bound for the multicast capacity 
is given by 
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where = | log (1 + x). For the same network, we can achieve all rates up to 



mm 

ser 



1 



X] «e Pu,v 

A(v) 



+ 1 



maxP^ 

As(v) 



(3) 



where = max{x,0}. Furthermore, the gap between the upper bound and the achievable rate is 
bounded by 

E log (1^(^)1)- (4) 

vev\{i} 

Remark 1: Note that, in the equal power case, i.e., Pu^y = P, the achievable multicast rate ([3]) has terms 
in the form of log (-^ + P) for some integer K > 1. Similar forms of achievable rate were observed in 
ifTOll , ifTSll , lfT6ll . |[25]l for some equal power Gaussian networks. 

The following subsections are devoted to proving Theorem [H 



A. Upper bound 

The cut-set bound |[27ll of the network is given by 



R < max min / {Xsy] Ysc\Xscy) . 

p{xv,v) 



(5) 



Though the cut-set bound is a general and convenient upper bound for the capacity, it is sometimes 
challenging to compute the exact cut-set bound in a closed form. This is due to the optimization by the 
joint probability density function (pdf) p{xvy)' In some cases, such as the finite-field networks in [[5]], 
ISl, [fT2ll . ifTSll . it is easy to compute the cut-set bound because a product distribution maximizes it. For 
the Gaussian case, however, the optimizing distribution for the cut-set bound is generally not a product 
distribution. 

Thus, we consider another upper bound which is easier to compute than the cut-set bound. This bound 
is referred to as the relaxed cut-set bound and given by 



R < min max / {Xsy] Ysc\Xscy) . 

SeT p[xv,v) 



(6) 



Due to the max-min inequality, the relaxed cut-set bound is looser than the original cut-set bound ((S]). 
For the relay network with interference, we can further simplify ^ as 



H^sy] Ys^\Xscy) 



= H^s^s^'iYs^lXscy) 
= H^s.s-i^sA^s-y), 

where the second and the third equalities follow by the structure of the network, i.e., 

• Xs^s {Xs^s^^, Xscy) Ysc, 



X 



s^y 



Yc 



. X 



X 



For cut S, the mutual information I{Xs^s^]Ysc\Xscy) is maximized when there is a perfect coherence 
between all inputs to a Gaussian MAC across the cut. Thus, we have 



m^^ I{Xs^s^]Ysc\Xscy) = V 

p(xv,V) 

Then by © and (|7]), the upper bound Q follows. 
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Fig. 2. Example: two-dimensional lattice constellation. 



B. Lattices and nested lattice codes 

Before proving the achievable part of Theorem [H let us establish some preliminaries for the lattices 
and nested lattice codes, which are key ingredients of our achievability proof. For a more comprehensive 
review on lattices and nested lattice codes, see lfT9ll , |[20ll , |[23l . An n-dimensional lattice A is defined as 
a discrete subgroup of Euclidean space with ordinary vector addition. This implies that for any lattice 
points A, G A, we have A + A^ G A, A — A^ G A, and G A. The nearest neighbor lattice quantizer 
associated with A is defined as 



Q(x) 



argmm ||x 



Al 



and the mod A operation is 



X mod A = X — (5(x) 



The (fundamental) Voronoi region of A, denoted by IZ, is defined as the set of points in closer to the 
origin than to any other lattice points, i.e., 

7^ = {x : g(x) = 0}, 

where ties are broken arbitrarily. In Fig. [2l an example of a two-dimensional lattice, and its Voronoi region 
are depicted. 

We now define some important parameters that characterize the lattice. The covering radius of the 
lattice r^"^^ is defined as the radius of a sphere circumscribing around IZ, i.e., 

r'^^ = min {r : 7^ C ri3}, 

where B is an n-dimensional unit sphere centered at the origin, and, thus, rB is a sphere of radius r. In 
addition, the effective radius of A, denoted by r^^, is the radius of a sphere with the same volume as TZ, 
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where Vol(-) denotes the volume of a region. The second moment per dimension of A is defined as the 
second moment per dimension associated with IZ, which is given by 



a\n) = -4:— ■ - / llxll^dx. 



In the rest of this paper, we also use Vol(A) and (t^(A), which have the same meaning as Vol(7^) and 
cr^(7^), respectively. Finally, we define the normalized second moment of A as 

(Vol(7^))'/" 

For any A, G(A) is greater than which is the normahzed second moment of a sphere whose dimension 
tends to infinity. 

Goodness of lattices 

We consider a sequence of lattices A^. The sequence of lattices is said to be Rogers-good if 

^cov 

lim — g = 1, 

which implies that A^ is asymptotically efficient for sphere covering |[20ll . This also implies the goodness 
of A^ for mean-square error quantization, i.e., 

lim G(A^) = — . 

n^oo 2716 

We now define the goodness of lattices related to the channel coding for the additive white Gaussian noise 
(AWGN) channel. A sequence of lattices is said to be Poltyrev-good if, for Z ^ A/'(0, a^I), 

Pr{Z ^ 7^} < e-^^^(^\ (8) 

where Ep{') is the Poltyrev exponent [22J and ji is the volume-to-noise ratio (VNR) defined as 

_ (Vol(7^))^/" 

Note that ([8]) upper bounds the error probability of the nearest lattice point decoding (or equivalently, 
Euclidean lattice decoding) when we use lattice points as codewords for the AWGN channel. Since 
Ep{ii) > for /X > 1, a necessary condition for reliable decoding is /i > 1. 

Nested lattices codes 

Now we consider two lattices A and A^. Assume that A is coarse compared to A^ in the sense that 
Vol(A) > Vol(Ac). We say that the coarse lattice A is a sublattice of the fine lattice A^ if A C Ac and 
call the quotient group (equivalently, the set of cosets of A relative to A^) A^/A a lattice partition. For 
the lattice partition, the set of coset leaders is defined as 

C = {Ac mod A} ^ {Ac n 7^}, 

and the partitioning ratio is 

: ^Vol(A) 



VVol(Ac), 

Formally, a lattice code is defined as an intersection of a lattice (possibly translated) and a bounding 
(shaping) region, which is sometimes a sphere. A nested lattice code is a special class of lattice codes, 
whose bounding region is the Voronoi region of a sublattice. That is, the nested lattice code is defined in 
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terms of lattice partition A^/A, in which A^ is used as codewords and A is used for shaping. The coding 
rate of the nested lattice code is given by 

-log|C| = log p. 

n 

Nested lattice codes have been studied in many previous articles [T8], [T9l, l[23]l . O, and proved to have 
many useful properties, such as achieving the capacity of the AWGN channel. In the next subsection, we 
deal with the nested lattice codes for the achievability proof of Theorem [B 

C. Nested lattice codes for a Gaussian MAC 

As an achievable scheme, we use a lattice coding scheme. In [TO], lattices were also used to prove an 
achievable rate of Gaussian relay networks with interference (called Gaussian MAC networks). However, 
they used the lattice as a source code with a distortion and then related the achievable distortion to the 
information flow through the network. Our approach is different from [lOJ in that we use lattices to produce 
coding and shaping lattices, and form nested lattice codes. As a result, our approach can handle unequal 
power constraints where incoming links have different power at a MAC. Our scheme is a generalization 
of the nested lattice codes used for the Gaussian two-way relay channel in [il5il . [1I61I . 

Let us consider a standard model of a Gaussian MAC with K input nodes: 

K 

y = 5]x, + z, (9) 

where Z denotes the AWGN process with zero mean and unit variance. Each channel input Xi is subject 
to the average power constraint P^, i.e., ^ Yl^=i{^f^Y ^ Pi- Without loss of generality, we assume that 
A > A > • • • > Px. 

The standard MAC in ([9]) is a representative of MACs in the Gaussian relay network with interference. 
Now, we introduce encoding and decoding schemes for the standard MAC. Let us first consider the 
following theorem which is a key for our code construction. 

Theorem 2: For any Pi > P2 ^ • • • > > and 7 > 0, a sequence of n-dimensional lattice chains 
A^ C A2 C • • • C A^ C Ag< exists that satisfies the following properties. 

a) A^, I <i < K, are simultaneously Rogers-good and Poltyrev-good while Ag< is Poltyrev-good. 

b) For any 5 > 0, - 5 < ^^(A^) < Pi, I <i < K, for sufficiently large n, 

c) The coding rate of the nested lattice code associated with the lattice partition Ag</A^ can approach 
any value as n tends to infinity, i.e., 

Rk = - log \Ck\ = 7 + On{l), 

n 

where Ck = {Ag< mod A^} and o^(l) ^ as n ^ oc. Furthermore, for 1 < z < X — 1, the coding rate 
of the nested lattice code associated with Ag</A^ is given by 

R,^- log ICI = Rk + ]: log ( ^ ) + 0,(1), 
n 2 \Pk J 

where Ci = {AJ mod A^}. 

Proof: See Appendix |Al ■ 
A conceptual representation of the lattice chain and the corresponding sets of coset leaders are given 
in Fig. [3] for a two-dimensional case. 

Encoding 

We consider a lattice chain as described in Theorem [2l We assign the i-th input node to the MAC with 
the set of coset leaders C^. For each input node, the message set {l, . . . , 2^^^} is arbitrarily mapped onto 
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Fig. 3. Example of lattice chain and sets of coset leaders. 



Ci. We also define random dither vectors ^ \]mi{TZi), I < i < K, where TZi denotes the Voronoi 
region of (we dropped the superscript for simphcity). These dither vectors are independent of each 
other and also independent of the message of each node and the noise. We assume that each is known 
to both the i-th input node and the receiver. To transmit a message that is uniform over {l, . . . , 2^^*}, 
node i chooses G Ci associated with the message and sends 

X, = (W, + U,) mod A,. 

Let us introduce a useful lemma, which is known as the crypto-lemma and frequently used in the rest 
of this paper. The lemma is given in [|23ll . and we repeat it here for completeness. 

Lemma 1 (Crypto-lemma [23]): Let C be a finite or compact group with group operation +. For 
independent random variables a and b over C, let c = a + 6. If a is uniform over C, then c is independent 
of b and uniform over C. 

By Lemma [B X^ is uniformly distributed over TZi and independent of W^. Thus, regardless of W^, the 
average transmit power of node i is equal to cr^(A^), which approaches Pi as n tends to infinity. Thus, 
the power constraint is met. 

Decoding 

Upon receiving Y = Xljli + where Z is a vector of i.i.d. Gaussian noise with zero mean and 
unit variance, the receiver computes 



K 



aY - ^ U^- mod Ai 



K 



K 



5](W,+U,)modA,-5]X, 



K 



K 



+ a ^ X^- + aZ - ^ U^- 
= (t + Z) modAi, 



mod Ai 
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where 



T = 



K 



5](w,-g,(w, + u,)) 



mod Ai 



K 



Wi + J](w,-g,-(w, + u,-)) 



mod Ai, 



(10) 



K 



Z = -(l-a)5]X, 



< a < 1 is a scaling factor, and Qj{') denotes the nearest neighbor lattice quantizer associated with 
Aj. We choose a as the minimum mean-square error (MMSE) coefficient to minimize the variance of the 
effective noise Z. Thus, 



a 



p. 



and the resulting noise variance satisfies 

1 



E 



n 



< 



p. 



(11) 



Note that, though the relation in ([TT]) is given by an inequality, it becomes tight as n ^ oc by Theorem 
[2l By the chain relation of the lattices in Theorem [2l it is easy to show that T G Ci. Regarding T, we 
have the following lemma. 

Lemma 2: T is uniform over Ci and independent of Z. 

Proof: Define W = Yf^=2^^ j - Qj^^ j + ^ j)) modAi, and, thus, T = (Wi + w) modAi. 

Note that W is correlated with X^, 2 < i < and Z. Since Wi is uniform over Ci and independent 
of W, T is independent of W and uniformly distributed over Ci (crypto-lemma). Hence, if T and Z are 
correlated, it is only through Wi. However, Wi and Z are independent of each other, and, consequently, 
T is also independent of Z. ■ 
The receiver tries to retrieve T from Y instead of recovering W^, I < i < K, separately. For the 
decoding method, we consider Euclidean lattice decoding [fT9ll -[231. which finds the closest point to Y 
in Ac. From the symmetry of the lattice structure and the independence between T and Z (Lemma©, 
the probability of decoding error is given by 



= Pr jZmodAi ^ Uc 



(12) 



where Qc{-) denotes the nearest neighbor lattice quantizer associated with and IZc denotes the Voronoi 
region of Ac- Then, we have the following theorem. 
Theorem 3: Let 



1 



log 



Pi 



y-^ p. 



+ Pi 



For any Ri < R^ and a lattice chain as described in Theorem |2l with Ri approaching Ri, i.e., Ri 
^1 + On (1), the error probability under Euclidean lattice decoding (fT2)) is bounded by 



Pe < e 



-n(£;p(22(«I-«l)] 
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Proof: See Appendix [Bl 



According to Theorem [3l the error probabihty vanishes as n ^ oc if < i?^ since > for 

X > 1. This imphes that the nested lattice code can achieve any rate below R[. Thus, by c) of Theorem 
Hand Theorem [3l the coding rate Ru I <i < K, can approach i?* arbitrarily closely while keeping 
arbitrarily small for sufficiently large n, where 



Remark 2: In theorem [3l we showed the error exponent of lattice decoding and the achievability of i?i 
directly followed. However, if we are only interested in finding the achievability of not in the error 
exponent, we can use the argument on the bounding behavior of lattice decoding in [l2T]l . which gives the 
same result in a much simpler way. 

Remark 3: Since Pi > • • • > Px, we have i?i > • • • > i?^- Now, consider the case that, for some 
I < K, the rates i?*, i + 1 < i < are zero while i?*, 1 < i < are nonzero. In this situation, nodes 
i + 1, . . . , K cannot transmit any useful information to the receiver, and, thus, we can turn them off so as 
not to hinder the transmissions of nodes 1, . . . , i Then, the variance of Z decreases and we have extended 
rates given by 



However, for the ease of exposition, we do not consider the transmitter turning-off technique and assume 
that nodes i + 1, . . . , K just transmit = when their coding rates are zero. 

D. Achievable multicast rate 

We consider B blocks of transmissions from the source to destinations. Each block consists of n channel 
uses. In block G {1, . . . , B}, din independent and uniform message W[k] G {1, . . . , 2^^} is sent from the 
source node s. It takes at most L = B+\V\—2 blocks for all the B messages to be received by destination 
nodes. After receiving L blocks, destination nodes decode the source message W = (Vl^[l], . . . ,W[B]). 
Thus, the overall rate is ^i?, which can be arbitrarily close to R by choosing B sufficiently large. 

Time-expanded network 

For ease of analysis, we consider the B blocks of transmissions over the time-expanded network |[2l], 
[El, ^TE, obtained by unfolding the original network Q over L + 1 time stages. In ^te, node ^ G 1/ at 
block k appears as v[k], and v[k] and v[k'] are treated as different nodes if k ^ k' . There are a virtual 
source and destination nodes, denoted by 5te and rfxE, respectively. We assume that 5te and s[k\^ are 
connected through virtual error-free infinite-capacity links, and, similarly, ^te and d[k\^ are. For instance, 
the network in Fig. [T] is expanded to the network in Fig. |4l Dealing with the time-expanded network does 
not impose any constraints on the network. Any scheme for the original network can be interpreted to 
a scheme for the time-expanded network and vice-versa. In our case, the transmissions of B messages 
W[k], k = 1, . . . ,B, from s io d ^ D over Q correspond to the transmission of a single message W from 
5te to (ixE ^ P^TE over ^te, where L^te denotes the set of virtual destination nodes. 

A main characteristic of the time-expanded network is that it is always a layered network [5J which 
has equal length paths from the source to each destinatiorQ. We define the set of nodes at length k from 
the virtual source node as 




(13) 




^Another characteristic is that the time-expanded network is always acyclic \2\. 
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Fig. 4. Time-expansion of the network in Fig. [T] 



and call it the k-ih layer. We use the subscript 'te' to differentiate parameters of Q and Qte- The set of 
nodes and edges of Qte are defined as 



^TE = {^Te} U i^TE U (^^uVte W j , 

^TE = {{u[k],v[k + l]):{u,v)eE,k = l,..., L} 
U{{s[k-l],s[k]) :k=l,...,L} 
U{{d[k],d[k+1]) : fc= 1,...,L}, 

where we define ^[0] = ste and d[L + 2] = rfxE- Note that, since Qte is layered, edges only appear 
between adjacent layers. From Vte and E^e, the other parameters, e.g., Ate(-)' 6)te(-)' ^te, Ste^ ^te, 
and Ate,s'(-)' similarly defined as A(-), B(-), S, S, F, and ^s{')^ respectively. 

Encoding 

We apply the nested lattice codes in Section IIII-CI over the all Gaussian MACs in the network. Thus, 
node v[k] is assigned with sets of coset leaders Cy[k],w[k+]^ w[k^] G QTE{y[k]), where k^ = k + 1. We do 
not change the lattice scheme over blocks, and, thus, C^[/e],^[/c+] = Cy^^j 

At node s[k], the indices {l, . . . , 2^^} are uniformly randomly mapped onto vectors in C^^^, w G ©(5). 
We define the random mapping as fs[k],w[f^+]{')' Then, node s[k] receives W = (VF[1], . . . , VF[5]) from 
s[k~] through the error-free link, where k~ = k — 1, and transmits 

on channel {s[k],w[k^]) using a random dither vector Us[k],w[k^]' ^odt v[k] that is not s[k] or d[k], 
the received signal is given by 

Yy[k] = [Tyik] + Zy[k]) mod A^, (14) 
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where 



T^v[k] — 



[k-],v[k] 



Qu,v (^^^ u[k~],v[k] + Vu[k-],v[k])) 



u[k-]e 

,ATE(^[fc]) 



mod A^, 



(15) 



and Zy[k] is an effective noise vector. In ([T4h . denotes the lattice associated with the incoming channel 
to node v with the largest power. Then, Ty[k] is decoded using Euclidean lattice decoding, which yields 
an estimate Ty^k]- Next, Ty[k] is uniformly and randomly mapped onto vectors in C^,^, w G Q{v). This 
mapping is denoted by fv[k],w[k+]{')^ and node v[k] transmits 



^v[k],w[k^] ~ fv[k],w[k^] \^v[k] 

on channel {v[k],w[k^]) using a random dither vector Uv[k],w[k+]' Node d[k], d e D, receives Y^ik] and 
computes T^ik]- It also receives (T^ii], . . . ,Td[k-]) from d[k~] through the virtual error-free infinite- 



G Ete are done independently. 



capacity link and passes ( T^[i], . . . , Td[k] ) to node d[k^]. 



We assume that all the random mappings fu[k],v[k+]^ 
Decoding 

While decoding, a virtual destination node rfxE ^ ^te assumes that there is no error in decoding 
Ty[k]'s in the network and that the network is deterministic. Therefore, with knowledge of all deterministic 
relation^ ([T5]) in the network, node ^te decodes W by simulating all 2^^^ messages and finding one 
that yields the received signal T^^j, = (Td[i], • • • , T^d[L+i]j • 

Calculation of the probability of error 

In the above decoding rule, we will declare an error if at least one of the following events occurs. 

• £i\ there is an error in decoding Tv[k] at at least one node in the network. 

• ^2- a message W ^ W exists that yields the same received signal T^^^, which is obtained under 
W, at at least one virtual destination node rfxE ^ ^te- 

Thus, the error probability is given by 

Pe = Pr{<?i U 82} 

<Pr{^i} + Pr{^2|fr}. (16) 
Let us consider the first term in ([T6l) . Using the union bound, we have 

Pr{fi} < 5Z 5] P^Ak]^ 



k=2 v[k]ev[k] 

\{s[k]} 



where 

Note that the summation is from k = 2 since nodes in the first layer do not have any received signal 
except for node ^[1]. By Theorem [3l at node v G 1^ \ {1} for any e > 0, Pe,v[k] is less than ^^^^ for 
sufficiently large n if 

Ru,v ~~ log iCii^x;! 



n 



+ 1 • Pu^v 



(17) 



^It is assumed that the all random dither vectors are known to destination nodes. Thus, ([Tsl) is deterministic. 
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for all u G A(^). Therefore, in this case 



Now, we consider the second term in (1161) . Under the condition ff, we have Tv[k\ = T^v[k\^ and, thus, 
the network is deterministic. Let us use the notation ^ u[k-],v[k\{W) and to explicitly denote the 

signals under message W. We say that node v[k] can distinguish W and if T^[/e](VI^) 7^ Ty[k\{W^)' 
Thus, from the argument of a deterministic network in [5J, the error probability is bounded by 

^r{8,\£l} < 2"^« ■ Pr I U {T,,,(W^) = T,,,(W^')} 

k ^TE 

= 2"^^^ • ^ PrjNodes in 5'te can distinguish W, W' , and nodes in 5'te cannot}. (18) 

-^TE^ 

Pte 

We briefly denote the probabilities in the summation in (fTSl) as 

Pr{P = 5TE,^ = 5^E}. 

Here, we redefine the cut in the time-expanded network ^te for convenience sake. From the encoding 
scheme, since the source message propagates through nodes s[k], k = 1, . . . ,L + 1, they can clearly 
distinguish W and W\ Similarly, if a virtual destination node rfxE cannot distinguish W and W\ nodes 
d[k], = 1, . . . , L + 1 cannot either. Thus, when we analyze the error probability (fTSl) . we can always 
assume that s[k] G Ste and d[k] G S'^g, /c = l,...,L + l, without loss of generality. 

From the fact that Qte is layered, we have 

Pr {V = Ste.V = =Pt{V = Ste.V = S'^^[l]} 



1[Ft{V = S^^^[k]\V = STE[k-lV = S^^^[k-]} 

k=2 

<1[Pt{V = S'^^[k]\V = STEnV = S'^^[k-]} , (19) 

k=2 

where S'TE[fc] and -STEffc] denote the sets of nodes in Ste and S^^ in the fc-th layer, i.e., 

STE[k]=STEnVTE[k], 

S^E[k] = S'TEr^VTE[k]. 
Also, from the fact that the random mapping for each channel is independent, we have 

Pr {f> = S'T,^[k]\V = STE[k-],V = S'^^[k-]} 
= H Pv {f> = {v[k]}\V = STE[k-],V = S^j.^[k-]} . (20) 

v[k]e 

Then, we have the following lemma. 

Lemma 3: Consider the time-expanded network Qte with independent uniform random mapping at 
each node. For any culEI Ste in Gte^ we have 



Pr {V = {v[k]}\V = STEik'iV = S^TEik']} < 2 
'From the definition, s[k] € Ste and d[k] e S^e^ fc = 1, . . . , L + 1. 



max Ru,v 
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for node v[k\ e S^^^lk], where S'^eW - -^te ^ ^TEik]- For node v[k] e S'^gfA;] \ ^^e[^]' 

we have 

Pr {V = {v[k]}\V = STE[k-],V = S^j^^ik-]} = 1. 



Proof: See Appendix O 
Thus, by ([IB-dini) and Lemma S it follows that 



Pr{£2\£t} < 2"^^ • IFte 



L + 1 

-n min J] E 

-Ste^ fc=2 y[k]e 



max Ru,v 
u[k-]e 



^ 2 ^TE Sf^^[k] V^TE,S(^['^]) J 



We now consider the following lemma. 

Lemma 4: In the time-expanded ^te with L + 1 layers, the term in the exponent of ([2T1) 



EE 



max A, 



is upper bounded by 



and lower bounded by 



mm 

k=2 _v[k]e yAxE,s(^[fe]) 

Srpg [k] 



L • min > ( max Ry ^ ] , 
(L — Irl + 2) • min > I max ^ \ . 



Proof: See Appendix [Dl 
Therefore, by (fTTl) . ([2T1) . and Lemma |4l Pr{52|^i} ^^^^ ^han | for sufficiently large n if 



L-|r| + 2 . ^ 

it < • mm > 

B s^v ^ 



1 



Thus, the total probability of error (fT6|) is less than e, and the achievability follows from (|22)) . 



(21) 



(22) 



L. G<2/7 between the upper and lower hounds 

To compute the gap between the upper bound Q and the achievable rate ([3]), we can rely on the 
following lemmas. 

Lemma 5: Assume that Pi > • • • > > 0. For any nonempty set ^4 C {1, . . . , X} and / = min74, 
we have 



Lemma 6: 



- log 
2 ^ 



y-^ p. 



+ l\Pi 



< log K. 



min{ai, . . . , a/,} - min{6i, ...,bk} 
< max{(ai - bi), . . . , (a^ - bk)}. 



The proof of Lemma |5] is given in Appendix El and the proof of Lemma |6] is omitted since it is 
straightforward. Using Lemmas |5] and |6l the gap in & directly follows. 
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Fig. 5. Linear finite-field symmetric MAC. 



IV. Linear finite-field symmetric networks with interference 

Let us consider a particular class of discrete memoryless relay networks with interference. The linear 
finite-field symmetric networks with interference are characterized by a special structure of MACs in 
them, which is shown in Fig. \5\ In more detail, the linear finite-field symmetric network with interference 
is described as follows: 

• Every input alphabet to a MAC at node v is the finite field, F^. 

• The received symbol at node v, yj^^ , is determined to be the output of a symmetric discrete memoryless 
channel (DMC) {¥q,p{yy\xy),yy) with input 

ueA{v) 

where Pu^y ^^q\ {0} denotes the channel coefficient. For the definition of the symmetric DMC, see 
[ESI Sec. 4.5]. 

• The input field size q and channel transition function p{yv\xy) associated with node v need not be 
identical. 

A major characteristic of the symmetric DMC is that linear codes can achieve the capacity [[281 Sec. 
6.2]. Using this, Nazer and Gastpar llllll showed that the computation capacity for any linear function of 
sources can be achieved in the linear finite-field symmetric MAC in Fig. [51 Also, in [T^, [fT3l], it was 
shown that linear codes achieve the multicast capacity of linear finite-field additive noise and erasure 
networks with interference, which are special cases of the class of networks stated above. Extending this 
line, we characterize the multicast capacity of the linear finite-field symmetric network with interference. 

Theorem 4: The multicast capacity of a linear finite-field symmetric network with interference is given 

by 

min > Cy, 

where Cy is the capacity of the channel (Fg,p(^^|x^), 3^^). 

The proof of Theorem [4] is very similar to the proof of Theorem [H The difference is that we use linear 
codes instead of the nested lattice codes. We show the outline of the proof in the next subsections. 

Remark 4: The capacity proof for linear finite-field additive noise networks in [12.1 can also be extended 
to the linear finite-field symmetric networks in Theorem [4l However, the proof in [fT2l relies on algebraic 
network coding, and, thus, it has a restriction on the field size, i.e., q> \D\.\n our proof, we do not use 
the algebraic network coding, and the field size is not restricted. 
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A. Upper bound 

As in the Gaussian case in Section Ull-Ai the upper bound follows from the relaxed cut-set bound 
In particular, for the linear finite-field symmetric network with interference, we have the Markov chain 
relation (X^^^c, Xscy) Xgc Ysc, where Xgc = {Xy : v G S^}. Using the data processing inequality, 
we have 

< I{Xsc] Ysc). 

Thus the upper bound is given by 

R < min max /(X^^^c; Ysc\Xscy) 
Ser p[xv,v) 

< min max I{Xsc] Y^c) 



in 



mm 

veS' 



B. Achievability 

Let us denote the vectors of channel input and output of the symmetric DMC (Fg,p(^^|x^), 3^^) as 



x.(l) 

y^v , • • • , 



T 



and 



J. , . . . , J. 



T 



, respectively. Without loss of generality, we assume 



that the encoder input is given by a uniform random vector G Fg^^''^ for some i?^ < 1. Then we 
have the following lemma related to linear coding for the DMC. 

Lemma 7 (Lemma 3 of [11]): For the symmetric DMC (Fg,p(^^|x^), 3^^), a sequence of matrices G 
^nx\nR^\ associated decoding function gy{') exist such that when = F^W^, Vi{g{Yy) ^ W^} < e 
for any e > and n large enough if = R'^ logg < C^. 

We now consider linear encoding for nodes in the network. We let 

X = d~^Y W 

and thus, 

X-y ^ ^ Pu,v'^u,v F-yT-y, 

where 

By Lemma Ul a linear code with sufficiently large dimension exists such that node v can recover with 
an arbitrarily small error probability if Ry < Cy. Now, we can do the same as in Section UlI-DI with (123)) 
replacing ([T5]) . and the achievability part follows. 

V. Conclusion 

In this paper, we considered the multicast problem for relay networks with interference and examined 
roles of some structured codes for the networks. Initially, we showed that nested lattice codes can achieve 
the multicast capacity of Gaussian relay networks with interference within a constant gap determined by 
the network topology. We also showed that linear codes achieve the multicast capacity of linear finite- 
field symmetric networks with interference. Finally, we should note that this work is an intermediate step 
toward more general networks. As an extension to multiple source networks, we showed that the same 
lattice coding scheme considered in this work can achieve the capacity of the Gaussian two-way relay 
channel within | bit [fTS]] . ifTTll . As another direction of extension, we can consider applying structured 
codes to networks with non-orthogonal broadcast channels. There is a recent work on the interference 
channel [26J which is related to this issue. 
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Appendix 

A. Proof of Theorem |2] 

Consider a lattice (more precisely, a sequence of lattices) with cr^(A^) = Pi, which is simultaneously 
Rogers-good and Poltyrev-good (simultaneously good shortly). In f20], it was shown that such a lattice 
always exists. Then, by the argument in [24J, we can find a fine lattice A2 such that A^ C and A2 is 
also simultaneously good. We let the partitioning ratio be 

/ Vol(A£) \ - ^ ^ / 1 \ ^ 

VVol(A^); \P2-S') \27reG{A^)J ^ ^ 

for some 5' > 0. Since the partitioning ratio can approach an arbitrary value as n tends to infinity, for 
any S > 0, n' exists such that we can choose S' < 5 when n > n'. We now have 

a2(A^)=G(A^).Vol(A^)t 

= G(A^) . 27re(P2 - S'), 

where the second equality follows from ([24h . Since is Rogers-good, n^^ exists such that 1 < 27reG(A2 ) < 
p^^^, , for n > Thus, for n > max{n^, n^^}, we have 

P2-S<a\A^)<P2. 

By repeating the same procedure, we obtain a lattice chain A^ C A2 C • • • C A^, where A^, I < i < K, 
are simultaneously good and Pi — S < cr^(A^) < for sufficiently large n. 

Moreover, by Theorem 5 of [19J, if A^ is simultaneously good, a Poltyrev-good lattice Ag< exists such 
that A^ C Ag< and the coding rate Rk can be arbitrary as n ^ oc, i.e., 

«-^-(^)— 

Given the coding rates i?^, 1 < i < X — 1, are given by 



1 






lOi 


n 




1 






lOi 


n 




1 








2 




1 








2 





VVol(AJ) 
/Vol(An\ 
ivoi(AJ)j+^'^ 

^ ) + + On(l), 

P 



where the third equality follows by the fact that A^ and A^ are both Rogers-good, and the fourth follows 
by the fact that ^^(A^) = Pi - o^(l). □ 

B. Proof of Theorem \3\ 

Let r^^^ and rf^ denote the covering and effective radii of A^, respectively. Then the second moment 
per dimension of r^^^B is given by 

' ' n + 2 
Next, we define independent Gaussian random variables 

Zi~A^(0,(T2l), 1 = 1,. ..,K, 
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and 

K 

Z* = (l-a)^Z^- + aZ. 

Then, we have the following lemmas. 

Lemma 8: The variance of Z*, each element of Z*, is denoted by Var(Z*) and satisfies 



K 



Var(Z*) = (l-o^)'5^^| 



2 ' 



< max 



' ^cov \ ^ p. 



J \ rf ' 



Lemma 9: The pdf of Z, denoted by Pz(x) satisfies 

Pz(x) < e''^f=^'^ •Pz*(x), 

where 

The above two lemmas are slight modifications of Lemmas 6 and 11 in lfT9l . The proofs also follow 
immediately from fT9]. 

Now, we bound the error probability by 

= Pr jZmod Ai ^ 7^c} 

<Pr|z^7^c} 

< e^^Ef^i . |2* ^ 7^c} , (25) 

where (1251) follows from Lemma [9l Note that Z* is a vector of i.i.d. zero-mean Gaussian random variables, 
and the VNR of relative to Z* is given by 

(Vol(Ac))^/" 
27reVar(Z*) 

> '^°""''':t' -°-w 

27re- ' 



22«i 27reG(A0 VEf=i^i 



+ Pi -0,(1) (27) 



where ([26b follows from Lemma [8] and the fact that A^, 1 < i < X, are Rogers-good, ([27l) from the 
definition of G(Ai), and (128)) from the fact that Ai is Rogers-good and Ri = + o^(l). When we 
consider the Poltyrev exponent, we are only interested in the case that ji > 1. Thus, from the definition 
of R\ and (1281) . we can write 

^ = 22(^1-^1) - o„(l), 
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for Ri < Rl. Finally, from (l25l) and by the fact that is Poltyrev-good, we have 

^ -n(£;p(22(«!-«l))-o„(l)) 



□ 



C. Proof of Lemma \3\ 

For notational simplicity, we prove this lemma in the standard MAC in Section IIII-C[ We assume that 
the uniform random mapping is done at each input node of the standard MAC, as was done in the network. 
Let A and A" be nonempty partitions of {1, . . . , K}, i.e., Au A" = {1, . . . , K}, and An A" = ^. We 
assume that A implies the set of nodes that can distinguish W and W, and A'^ implies the set of nodes 
that cannot. For node i e A, Wi{W) and Wi{W') are uniform over Ci and independent of each other 
due to the uniform random mapping. However, for node i G we always have Wi{W) = Wi{W'). 
Thus, if A = 0, T{W) = T{W') always holds, i.e., 

Pr{T{W) = T{W')\V = A,V = A''} = 1. 

If A ^ 0, given V = A and f> = A", the event T{W) = T{W') is equivalent to f{W) = f{W'), where 

fiW)= ^(W,(W^)-g^.(W,-W + U,-)) modAi, 

JeA 

and T{W') is given accordingly. Now, let / = min^, then 
T'{W) = f{W) modA; 

^liW) + (w,(H^) - Qj{W,{W) + U,)) 

\{l} 

which follows from the fact that Ai C A/, and ,thus, (xmod Ai) mod A/ = xmod A/. Note that, due 
to the crypto-lemma and the uniform random mapping, T^(VI^) and T\W') are uniform over Ci and 
independent of each other. Therefore, 

Pt{T{W) =T{W')\V = A,V = A''} = Pr|T(T^) =f{W')\V = A^ 

< Ft{T\W) = T\W')\V = A} 
1 



mod A 



Thus, by changing notations properly to those of the network, we complete the proof. 



□ 



D. Proof of Lemma |?1 

In the time-expanded network, there are two types of cuts, steady cuts and wiggling cuts 0. The steady 
cut separates the nodes in different layers identically. That is, for a steady cut aSte, ^[k] G S'te for some 
k if and only if ^[1], . . . , ^[L + 1] G S'te- Let us denote the set of all steady cuts as Fte- Then, since 

Tte ^ Tte? 

L+l / \ L+1 / 

min max Ruv\ ^ niin max Ruv 

k=2 y[k]^ yATE^^(v[fc]) J k=2 y[k]e \^Arj.^^g(v[k]) 



L • min >^ ( max Ruv) - 
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We now prove the lower bound. For any two cuts Si and 5^2 in i.e., Si, S2 G F, define that 

^{Si, S2) = ^ i max Ru^y j , 

where i?^,^^ = if {u, v) ^ Then, we have the following lemma 

Lemma 10: Consider a sequence of non-identical cuts S'l, . . . , S'l^ G F and define Sl'^i = 5^1. For the 
sequence, we have 

L' 

k=l k=l 

where for A; = 1, ... , L\ 

The proof of Lemma [10] is tedious but straightforward. Similar lemmas were presented and proved in 
[[5l Lemma 6.4], (Si Lemma 2], and the proof of Lemma [TOl also follows similarly. 
Now, since S'^ G F, it follows that 

^(^L > min ^ ( max i?^,^ j 

5^ f mjx Ru,, J . (29) 



mm 

ves^ 



Also, since S^TEi^l's correspond to cuts in V, we can rewrite 

L+l ( \ L+1 



V V max RuA =^min V^(5TE[fc 



mm 



Since there are |F| = 21^'"^ different cuts, at least the first L— |F|+2 of the sequence S'xEil], • • • , 5'te[^+1] 
form loops, and, thus, by Lemma [TOl and ([29b , we have 

min (5TE[fc-], 5teW) > (L - |F| + 2) • min V maxi?,,. 



Proof of Lemma [5| 

We first consider the case that 1 G yl, and the case that \ ^ A afterward 
a) 1 G A 

In this case, / = 1, and the gap is 



□ 
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b) 1 ^ A 

Since 1 ^ A, \A\ < K — 1. Now, the gap is given by 



<-log{l + {K-lfPi) 
<llog(l + (X-lf) 



^log Pi 



p. 



□ 
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